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Masses of the charm and bottom quarks are important inputs to pre¬ 
cision calculations of Higgs boson observables, such as its partial widths 
and branching fractions. They constitute a major source of theory uncer¬ 
tainties that needs to be better understood and reduced in light of future 
high-precision measurements. Conventionally, Higgs boson observables 
are calculated in terms of rUc and m^, whose values are obtained by av¬ 
eraging over many extractions from low-energy data. This approach may 
ultimately be unsatisfactory, since me and mb as single numbers hide var¬ 
ious sources of uncertainties involved in their extractions some of which 
call for more careful estimations, and also hide correlations with addi¬ 
tional inputs such as Aiming at a more detailed understanding of 
the uncertainties from me and in precision Higgs boson analyses, we 
present a calculation of Higgs boson observables in terms of low-energy 
observables, which reveals concrete sources of uncertainties that challenge 
sub-percent-level calculations of Higgs boson partial widths. 
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1 Introduction 


Precision studies of the Higgs boson are crucial for understanding the origin of elec- 
troweak symmetry breaking and searching for new physics beyond the Standard Model 
(SM). Many well-motivated new physics scenarios predict percent-level deviations 
from the SM for the Higgs couplings [T]. Accordingly, it is hoped that data at a 
later stage of the LHC or future facilities will enable (sub)percent-level determina¬ 
tions of the Higgs boson branching fractions and the important partial widths 
{Q = c,b) |21 El in E] • However, this does not necessarily imply that percent-level new 
physics effects can be probed. In fact, the power of precision Higgs analyses is limited 
by theory uncertainties, which are dominated by the parametric uncertainties from 
the input parameters, especially the charm and bottom quark masses [6l[7ll8]El In 
terms of the scale-invariant masses in the MS scheme [i.e. solutions to mQ(p) = /i], 
the uncertainty propagation, according to Ref. [7], 


lOMeV • 


^ Amfc(mfe) 

^H^bb ~ lOMeV 


X 0.56%, 


( 1 ) 


indicates that this parametric uncertainty is at the level of a few percent at present, 
if the input quark masses are taken from the PDG particle listings uni. 


mc{mc) = 1.275(25) GeV, rribirnb) = 4.18(3) GeV. (2) 


Note also that propagates into the calculations of all branching fractions since 

bb is the dominant decay channel of the SM Higgs boson. In light of the projected 
experimental precision on the Higgs observables, it is thus highly desirable to have a 
detailed understanding of this dominant theory uncertainty, and hence the improve¬ 
ments needed. In particular, instead of treating rric and nib as single numbers to be 
read from the PDG, we will make an initial attempt to incorporate the extraction of 
rUc, nib from low-energy experiments into precision Higgs analyses. As a result, the 
vague notion of “uncertainties from quark masses” will be decomposed into concrete 
sources, some of which represent a serious challenge that calls for further investigation 
of the low-energy observables. Further details of this work can be found in Ref. m- 


2 A global approach to precision Higgs analyses 

A conventional approach to precision Higgs analyses is to regard nic and nib as in¬ 
put “observables”, namely to put them on the same footings as mz and nin- As 
such they are assigned “experimental” central values and error bars, as in Eq. (1^ . 
However, unlike mz and mn, quark masses are not well-dehned observables due to 

*The perturbative uncertainty for T is well below 1%, thanks to the N"^LO calculation [9]. 
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confinement. Instead, they are just parameters in the SM Lagrangian. Their values 
are extracted from the true observables whose theory predictions depend on them. In 
fact, the PDG quark masses quoted in Eq. ([2]) are obtained by averaging over many 
results in the literature, which are extracted from a variety of observables. Much 
information is hidden by this averaging procedure, making the conventional approach 
of using the numbers in Eq. (|2]) unsatisfactory in several aspects. First, the averaging 
unrealistically assumes no correlations among the various quark masses extractions, 
some of which use similar data and/or methods. Second, there are correlations be¬ 
tween niQ and because enters the quark masses extractions as an input, but 
such correlations are not retained in Eq. ([2]). mg and as are then treated as inde¬ 
pendent inputs when the Higgs observables are calculated, which is strictly speaking 
incorrect. Furthermore, the meaning of the error bars in Eq. ([2]) is obscure. They 
contain not only various experimental and theoretical uncertainties associated with 
many different observables, but also a self-described inflation of uncertainties by the 
PDG (see the “Quark masses” review in Ref. HD]). The latter is introduced to ac¬ 
count for possibly underestimated uncertainties in some of the reported quark masses 
in the literature. The problem of uncertainty underestimation was hrst noticed in 
Ref. [12], and its impact on precision Higgs analyses will be discussed below. 

These unsatisfactory aspects of the conventional approach can in principle be 
eliminated if we include only true observables {O^} in the fit) where the function 


= 




of({i}) - or- 


ij 




-1 


OfiiXk}) 


^expt 


(3) 


with appropriately determined uncertainties and correlations contained in the covari¬ 
ance matrix V, is minimized with respect to the inputs of the calculation {X^}. rric 
and mb are in the set {X^}, but not in the set {O^}. 

For such a £t to be useful, the £t observables {O^} should include those con¬ 
tributing to the extraction of and mb- Interestingly, while our starting point is 
precision Higgs analyses, the observables dominating the PDG average of mg are as¬ 
sociated with much lower energy scales E ~ mg <C mi^. They include, for example, 
low [121 HSl HU HSl [16] and high [T7| [HI [191 l2D] moments of e+e“ —)■ QQ inclusive 
cross section, dehned by 



ds 

gn+l 


Rq^s), 


where Rq = 


cr(e+e —)■ QQX) 
cr(e+e“ —)■ ’ 


(4) 


variants of these moments pil |D2] , and moments of lepton energy and hadron mass 
distributions of semileptonic B decay [231 [2T[ [25]. 
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As a result, we are led to the following picture: 
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The low-energy observables {0]°"} play an important role in precision Higgs analyses, 
because they are sensitive to the same set of inputs as the Higgs observables 
The role of low-energy observables is not obvious in the conventional approach, where 
a large amount of information from has been highly processed into just two 

numbers rric and rrib. As we strive for higher-precision calculations, this information 
should be resolved, and the low-energy observables should be more directly engaged. 
In this way it is in principle straightforward to include correlations among the ob¬ 
servables, and retain the correct dependence in each observable. We propose this 
global approach as a long-term goal for the precision Higgs analysis program, which 
will become more relevant as experimental precisions on improve over time, 

for both rigorous tests of the SM and hts to SM extensions. The observables set can 
also be enlarged to include precision electroweak observables (e.g. Z-pole observables, 
mw, and LEP2 data) to make the global approach even more powerful. 


3 Anatomy of theory uncertainties in r 


H^hh 


To assess the theory uncertainties in calculating without performing a full- 

fledged £t, in particular the uncertainties from our imprecise knowledge of rric 
and mfc, it is helpful to pick out two low-energy observables from the set {0]°'"}, 
and eliminate nic and mb in favor of them in the functions mb, 

We will choose dehned in Eq. (|1]), as the two low-energy observables, 

motivated by the simplicity of their calculations and the small quoted uncertainties 
for me and m^ extracted from them. These moments can be calculated by the method 
of relativistic QCD sum rules [26] (see e.g. Refs. [23 [2H] for reviews )li], which relate 
them to the vector current correlators: 

127r2 / 


M9 = 


ni 


where 


- q^,qu)'^^Q{q‘^) = -i 


ngtf) , 

(6) 

q^=0 


a:e*’'yo|rj,p)jt(o)|o). 

(7) 


^We note in passing that the sum rules approach has been recently recast by the lattice QCD 
community [221 IMl 131] • See Ref. [S] for its possible impact on future precision Higgs analyses. 
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is the electromagnetic current of the quark Q{= c, b). IIq can be calculated via an 
operator product expansion: 

{n,)( ‘ In- ln‘ 


{2mQ{Hrn)y 


i,a,b 


K 


K 


( 8 ) 

where Qq is the electric charge of the quark Q, and are functions of the number 
of active quark flavors uj (4 for Ai^ and 5 for A4^). The two terms come from 
perturbation theory and nonperturbative condensates, respectively. Low moments 
(small n) are preferred so that the perturbative part, which has been calculated up to 
O (af) [32], dominates. The nonperturbative piece is dominated by the gluon 

condensate contribution, which has been calculated to next-to-leading order [33] , and 
is nonnegligible only for the charm quark. 

The moments are functions of mg and ag, which are renormalized at fim and 
/Xq, respectively, in the MS scheme. Physical observables like Ai^ should not depend 
on the renormalization scales /im, Pq. But when they are calculated to hnite order 
in perturbation theory there is residual scale dependence, which is conventionally 
used to estimate the effects of unknown higher-order terms. It used to be a common 
practice to set fim = l^a, hut it is argued in Ref. [12] that the perturbative uncertainty 
is generally underestimated in this way. Keeping and fia separate, we can invert 
Eq. ([8]) to extract mg, as has been done in Refs. [121 CH], 


rricimc) 


rricimc) as{mz), Ail, 1^1, Ml 


c,np 


(9) 


mfe(mfe) 


mb{mb) as{mz),Ml,i4n,iy 


( 10 ) 


Here we have assumed the input is at the scale mz, from which as{^a) can be 
derived by the renormalization group (RG) equations [51]. Also, RG equations allow 
us to convert the extracted mg{^rn) to mg{mg). It is clear that the quark masses 
depend not only on the low-energy observables M^, but also on as{mz) and the 
renormalization scales in the calculation of M^. All this dependence is retained in 
Eqs. dH]) and flTU]) . and is eventually propagated into the calculated Higgs observables 
in We will focus on the partial widths T ^^he following. 

Neglecting correlations between rric and (which is justihed given the large 

uncertainties in MY^'^, see Ref. m), we have 


T 


H^cc 


= T 


H^cc 


rricimc), Yh 


= T 


H^cc 


H^bb 


j. = T 


H^bb 


= T 


H^bb 


{ot}.Mi,Y^,Yc..YH.Ml 


{0“},mb(mb),/i^ 


{Ot},M\,Yn,YaMy 


c,np 


( 11 ) 


( 12 ) 
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where {O™} = {mz-, Gp, a{rnz), rrit, tuh} is the input observables set 

familiar in precision electroweak analyses, /i^ ~ O {rriH) collectively denotes the 
renormalization scales chosen in the calculations of Th^cc and and should not 

be confused with ~ O (2mQ). Note that the as{mz) dependence has changed 

in the second equalities in Eqs. flTTll and flT^ to account for the correlations between 
as and mq reflected in Eqs. (jO]) and ffTOj) . 

Eqs. JIB and f|12p show the decomposition of the theory uncertainties in Vh^cc) 
^H^bb- particular, what people usually refer to as “uncertainties from rric and mi," 
are broken down into concrete sources of uncertainties, the most important ones being 
the parametric uncertainty from the measurements of the low-energy observables 
and the perturbative uncertainty reflected by the residual scale dependence. 

While the former is straightforward to quantify and interpret ji| the latter necessarily 
involves artificial prescriptions that may lead to bias in its estimation, as we will 
discuss below. We also note that the parametric uncertainty from asijnz) is affected 
by the dependence of the extracted rngirng) on asijnz), and is found to be smaller 
than the incorrect estimate neglecting this correlation; see Fig. [2] below. 

To visualize the perturbative uncertainty from the dependence on /im, /^a, we £x 
all the input observables at their experimental central values listed in Ref. [n] and 
set hh = ttT'H, and make contour plots for the calculated Th^ce, ^H^bb 
plane. These plots, shown in Fig.lH illustrate the propagation of the scale dependence 
from the low-energy observables to the Higgs boson partial widths, as the latter are 
seen to depend on the renormalization scales chosen when calculating the former. The 
necessity to go beyond jjLm = l^-a is clear since the diagonal does not capture all the 
scale dependence. To estimate the perturbative uncertainty, a common practice is to 
vary the renormalization scales within a factor of 2 around a characteristic scale of 
the process. But this does not directly apply to and /Tq, because the receive 
contributions from all \/s; see Eq. (jl]). Therefore, for illustration we will vary /i^ 
and fia independently within an adjustable range [/imin, A^max], and plot the estimated 
perturbative uncertainty as a function of /imin for a few choices of /imax in Fig. [2l The 
leading parametric uncertainties are shown in the same hgure for comparison. 

It is seen from Fig. [2] that the perturbative uncertainty from /Xm, /^a is very sen¬ 
sitive to the somewhat arbitrary choice of /Xmin; and can dominate the total theory 
uncertainty for T h^cc, ^n^bl lower renormalization scales are allowed in the calcu¬ 
lations of A4^. This represents a serious challenge for higher-precision Higgs boson 
partial widths calculations, and can be truly overcome only by calculating Ai^, or 
equivalently (^) flQiq^)\^ 2 ^Q fo higher order. Before this calculation is done, ideas 
on more enlightened prescriptions for the uncertainty estimation may be helpful. Two 

tNote, however, a complication due to the fact that the “experimental” uncertainty of 
contains a contribution from perturbative QCD input for y/s > 11.2 GeV where no data is available 
at present. This leads to a large “experimental” uncertainty in AlJ, and explains why the second 
moment \s preferred for extracting mi,. The situation is expected to improve in the future. 
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Figure 1: Contours of T^H^bl MeV in the plane. Unlabeled con¬ 

tours represent decreasing values toward the top-left corner in steps of 0.002, 0.005, 
respectively. The diagonal pm = Pa does not capture all the perturbative uncertainty. 

possibilities, the convergence test [I6] and optimal scale setting |35l IMl EZ], are dis¬ 
cussed in Ref. m, but both of them are still unsatisfactory at present. It is possible 
that the actual situation is better in a global fit incorporating more observables 
in {01°®}, but it remains to investigate other low-energy observables sensitive to rric, 
nib and see if their calculations are plagued by similar scale-setting ambiguities. 

In passing we briefly comment on the perturbative uncertainty from which 
is usually estimated by varying hh from mn/^ to 2mH', see e.g. Ref. [7]. It can 
be argued that this procedure is incomplete, since the calculations of F h^cc, 
also involve more than one renormalized parameters, and the renormalization scales 
chosen for mg and as need not be equal. However, we have found that in this case, 
keeping the scales separate does not signihcantly increase the estimated uncertainty. 
Thus, the estimates in the literature with only one scale are expected to be robust. 

4 Conclusions 

The success of the SM is based on its agreement with data for processes across all 
accessible energy scales. Charm and bottom quark masses provide a bridge connecting 
two distinct scales mq and m,H- At the next level of precision in testing the Higgs 
sector of the SM, theory calculations need to be improved to match the projected 
experimental precision. For a better understanding and a more consistent treatment 
of the theory uncertainties, it is desirable to resolve the information contained in 
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^Jmin/GeV Mj,i„/GeV 


Figure 2: Percent relative uncertainties in Th^cci ^H^bl functions of /Xmin from 
various sources: perturbative uncertainty with /i^ax = 4 GeV, /i^ax = GeV (red 
solid) or alternatively p^ax = 3,5 GeV, /i^ax = 13,17 GeV (red dashed), parametric 
uncertainties from or M.\ (orange), as{mz) (cyan solid), (blue, for Th^cc 

only) and mn (purple). The parametric uncertainty from as{mz) incorrectly calcu¬ 
lated assuming no correlation with mg (cyan dotted) is also shown for comparison. 
See Ref. m for more details. 

rric and rrib extracted from low-energy data. We have proposed a global approach 
involving low-energy observables and Higgs observables as a long-term goal for the 
precision program, and performed a hrst calculation that connects the two sets of 
observables. The analysis gives us a more detailed understanding of the “uncertainties 
from rric, mb ", and points to future directions for the precision Higgs analysis program. 
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